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Abstract 

We calculated the complete next-to-leading order(NLO) QCD corrections to the ttbb pro- 
duction process at a 77 collider in the standard mode^^M) . The calculation of the one- loop 
QCD correction includes the evaluations of the hexagon and pentagon amplitudes. We stud- 
ied the NLO QCD corrected total cross section, the distributions of transverse momenta of 
final top- and bottom-quark states, and the dependence of the cross section on renormaliza- 
tion scale fi. It shows that NLO QCD correction generally increases the LO cross section 
in our chosen parameter space, and the K-factor varies from 1.70 to 1.14 when colliding 
energy goes up from 400 GeV to 2 TeV. We find that the correction distinctly changes 
the distributions of transverse momenta of the final top- and bottom-quark states, and the 
NLO QCD correction obviously improves the independence of the cross section for process 
77 — *■ ttbb on the renormalization scale. 



PACS: 14.65. Ha, 14.80.Bn, 12.38.Bx 



I. Introduction 



Recently, the report of the top-quark mass measurements from the CDF and DO experiments 
at Fermilab, presents the preUminary world average mass of the top quark is mt = 172.5 it 
1.3{stat) lb 1.9{syst) GeV, which corresponds to a 20% precision improvement relative to the 
previous combination [I]. Due to the huge top-quark mass being much heavier than the VF-boson, 
before the top-quark hadronization it undergoes dominantly the weak decay via t W~^b[2\, that 
has two important consequences different with other quarks, the narrow resonance around the 
energy 2mt is absent and the perturbative QCD is reliable to study all the threshold region. We 
believe the large top mass value will open up new vistas of electroweak physics and possibilities 
for probing the effects beyond the SM. For example, since the Higgs mechanism in the SM and 
other extended models predicts that the strength of quark-Higgs Yukawa coupling is proportional 
to the quark mass[3], one can measure top-quark Yukawa couplings in high precision to probe 
the SM or discover new physics. There have been many works concerning the study of the 
top-quark physics at colliders [H El El El E] . All those indicate that the precise study of the top 
physics is accessible. In the proposed and planned experiments, such as at the CERN LHC and 
the ILC, more interests are focused on the developing and understanding top physics and top 
characteristics. 

The future e~^e~ linear collider(LC) not only can be designed to operate in e~^e~ collision 
mode, but also can be operated as a 77 collider. This is achieved by using Compton backscattered 
photons in the scattering of intense laser photons on the initial e^e~ beams. The top-quark 
pair production at photon-photon collider is also useful for top-quark physics study. It has been 
found that the 77 — > tt production rate with high 77 colliding energy, is much larger than 
that from the direct e~^e~ — > tt production due to the s-channel suppression of later process 

At the future International e~^e~ Linear Collider (ILC), the colliding center-of- mass-system 
(c.m.s.) energy can be raised up to 2 TeV. At TeV scale high colliding energy, the final states 
will be very complex and expected to be in multi-particles or jets with large production rates, 
and top and Higgs signatures naturally are included in these processes. For these processes with 
large cross sections, a leading order analysis is not adequate to make detailed predictions for their 



cross sections. In fact, the processes with multiple final particles are particularly interesting, 
since such processes often proceed through one or more resonances that subsequently decay, 
or they represent an irreducible background to such resonance processes. For example, the 
associated production 77 ttH is an important process in probing top-Higgs Yukawa coupling 
in high precision, and was already studied in Ref.|12j. But after Higgs boson and top-quark 
decay, the ttH associated production process at 77 collider would have the same final states 
(W^bW'bbb) as the process 77 — > ttbb . Therefore, it is important to emphasize that the ability 
to distinguish top, Higgs boson or other new particle signatures at linear colliders, crucially 
depends on the understanding of the corresponding backgrounds of the corresponding processes 
with multi-particle final states. In order to exert all the abilities of future colliders, precise 
theoretical predictions including higher order corrections to multi-particle production processes 
are necessary. 

In previous work, people used " double-pole approximation" (DPA) [TB] to handle the evalua- 
tions of complete one-loop strong and electroweak calculation for process including four particle 
final states. Recently, the methods for the calculation of scalar and the tensor 6-point integral 
functions were provided [TU [TB]. With the approach provided in Ref. |14j . the complete elec- 
troweak corrections to the e~^e~ — > 4/ processes, which are relevant for W-pair production, was 
presented by A. Denner, S. Dittmaier, M. Roth, L.H. Wieders, and the results were compared 
with those in Ref.|17j by using DPA method. 

In this paper we present the calculations of the cross sections of the process 77 ttbb 
at the leading-order(LO) {0{al^a'j.)) and QCD next-to-leading-order(NLO) {0{al^al)). The 
paper is organized as follow: The tree-level analytical calculation of the cross section for the 
process 77 — > ttbb is given in section II. In section III the analytical calculation of the NLO 
QCD corrections is presented. The numerical results and discussions are given in section IV. 
Finally, a short summary is given in section V. 




II. Analytical calculation of the cross section for the process 

77 tibb at the tree-\eve\{0 {al^,a1)) 

The process 77 — > ttbb can be induced via t— and u— channel at the tree-level. In this section 
we consider the tree-level(C'(ag^Q^)) contribution to the cross section for the process 77 ttbb 
. We denote the process as 

liPi) + 7(P2) ^ t{p3) + i[pi) + b{p5) + Hpfi), (2.1) 

where the four- momenta of incoming photons are denoted as pi and p2, and p^, p4, ps, pQ 
represent the four-momenta of the final particles. The FeynArts3.2 package|l8j is adopted 
to generate tree-level Feynman diagrams and convert them to corresponding amplitudes. We 
present the t-channel tree-level Feynman diagrams involving strong interaction for the process 
77 — > ttbb in Fig.l. The u-channel diagrams, which can be obtained by exchanging initial 
photons of the corresponding t-channel ones, are not drawn there. 

The tree-level amplitude M. of the process 77 — > ttbb at the O^aewOs) order, is then expressed 
as below. 

10 t 

^ = EE^f- (2-2) 

i=l j=u 

where the amplitudes (A^f \ ^ = 1; •••) 10) correspond to the diagrams in Fig{TI^l-10) respectively. 



Then we get the differential cross section for the process 77 — > ttbb at the tree-level as 

datree = ^\M\^d<^^, (2.3) 

where d^4 is the four-body phase space element. The summation is taken over the spins and 
colors of final states, and the bar over the summation in Eq. ()2.3p recalls averaging over initial 
spin states. The calculation of the amplitudes of tree-level diagrams for 77 — > ttbb process is 
implemented by using FormCalc4.1 package|19j. The integration program of the four-body phase 
space is based on the FormCalc4.1 package, and created by using the factorization expression of 
four-body phase space element [2D], 

d$4(Pl +P2,P3,P4,P5,P6) = 77-dQ^d^>3(Pl + P2,P3,P4,Q)d^2{Q,P!i,P6)- (2.4) 

ZTT 

where Q = + pQ and the actual integration of three-body phase space element d^3{pi + 
P2,P3,Pi,Q) is parameterized using 2to3.F program in FormCalc4.1 package[19j. and the com- 
plete integration over the four-body phase space is performed using Monte Carlo integrator 
Vegas. The two-body phase space d^2iQ,P5,P6) is given as 

dMQ,P5,P,) = .^\2'^^^^l'Z'' —d^5d{cose,), (2.5) 

{2iTy 8(y^ 

where the kinematical function X{x,y,z) is defined by 

A(x, y, z) = + y"^ + — 2xy — 2yz — 2zx. (2-6) 

In order to check our four-body phase space integration program, we calculated the LO cross 
section of e"^e~ — > ud^^U^ process by adopting fixed-width method and the same input param- 
eters as used in Ref. [17J. In Table [1] we list the the LO cross sections obtained from Ref.|17j. 
by using GRACE2.2.0 system|21j, and by adopting our 4-body phase space integration pro- 
gram together with FeynArts3.2 and FormCalc4.1, respectively. It demonstrates the numerical 
integration results by using our created program are in good agreement with others. 

III. The NLO QCD corrections to the process 77 ttbb 

The FeynArts3.2 package is used to generate QCD one-loop Feynman diagrams of 77 — > ttbb 
process at the order of 0{aewOii), and then to convert them to corresponding amplitudes. The 





aLo{fb){Rcf.m 


aLo{fb){GRACE) 


o-Lo(/^) (ours) 


200 
500 


661.3(3) 
260.9(1) 


661.24(3) 
260.85(2) 


661.26(4) 
260.87(2) 



Table 1: The comparison of the numerical results of the LO cross section for e^e~ — > 
udiJ,~u^ process by using GRACE2.2.0 system and our integration program with those 
presented in Ref.[T7j. 

QCD one-loop Feynman diagrams can be classified into 140 self-energy diagrams, 172 triangle 
diagrams, 108 box diagrams, 48 pentagon diagrams and 12 hexagon diagrams. We also use the 
FormCalc4.1 package [TSTj to calculate the amplitudes of one- loop Feynman diagrams. But the 
original FormCalc4.1 package doesn't possess the function to calculate the amplitudes including 
6-point integrals, we have to create some codes in FormCalc4.1 in order to handle amplitudes 
relevant to hexagon diagrams. As a representative selection, we present the hexagon Feynman 
diagrams of the 77 ttbb process in Figl2j There exist both ultraviolet (UV) divergency and soft 
infrared(IR) singularity in the contribution part of virtual gluon one-loop diagrams for 77 — > ttbb 
process, but no collinear IR singularity due to the nonzero masses of top and bottom quark. 
Dimensional regularization(DR) scheme in D = 4 — 2e dimensions is adopted to isolate both 
IR and UV singularities. By adopting MS-scheme to renormalize the strong coupling strength 
and the OS'-scheme to renormalize the masses and fields of top- and bottom-quark, the UV 
singularities are vanished. 

In the total cross section, the soft IR divergency contributed by the virtual correction part 
should be cancelled by the contribution from real gluon emission process 77 ttbb +g at tvee- 
level{0 (a^^a^)) , which presents the same order contribution as the virtual correction does. The 
real gluon emission process is denoted as 

liPi) + 7(P2) ^ t{p3) + i{pi) + b{p5) + Kpg) + gipj). (3.1) 

We adopt the two cutoff phase-space slicing(TCPSS) method [22] to calculate the real gluon 
emission process. Since there is no collinear IR singularity, we introduce only an arbitrary small 
soft cutoff 5s to separate the 77 ttbb +g phase space into two regions, according to whether 
the energy of the emitted gluon is soft(£'7 < 5s\^s/2) or hard(i?7 > 5sy/s/2). After a lengthy 
calculation similar with those shown in references [231 [24t [25| [26] , we can get the expression of 
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Figure 2: The hexagon Feynman diagrams for 77 — > process. 



Cso/t for 77 — > ^i^ft process as, 



soft 



Pis 

2-K 



17 1 

o (f35 + 546) + ^ (536 + 545) - 7 (534 + 556) 




where (i,j = 3,4,5,6) are defined as, 



9^o{P^,Po) = (1^) r(l 



+ e) 
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i ln2((T,) + ^ ln"(cjj) + 2Li2 (1 - at) + 2Li2 (1 - aj) 
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(3.3) 



In above equation, X{sij,m^ ,mj) is the kinematical function defined in Eq. (|2.6p . AE = 6s\/s/2^ 



{p, + p,)\ m^=mi = mt, = = m, and 



Pi 



Sij + mf - m| 
1 - Pi 



(3.4) 



l + pi 

Our calculation shows in the total cross section the soft IR singularity induced by the one-loop 
virtual gluon correction is exactly cancelled by the IR singularity part offered by the soft gluon 
emission process 77 tibb (g). The hard gluon emission cross section a hard for -5-7 > (^s\/s/2 
is finite and can be calculated numerically in four dimensions by using Monte Carlo method. 



The integrations of 4-body phase space in evaluations of crtree, (^virtual and ago ft are carried 
out by using our created 4-body phase space integration program as described in last section. 
But in evaluation of ahard of 77 tibb +g process, we apply CompHEP-4.4p3 program[29j to 
calculate the tree-level amplitude and the integration of 5-body phase space. Finally, the NLO 
QCD corrected total cross sections for process 77 — > ttbb can be obtained as 

(^NLO = entree + (^virtual + CTsoft + CThard- (3.5) 

The QCD NLO corrected cross section gnlo is both UV- and IR- finite. All the UV and IR 
divergences are cancelled analytically. 

IV. Numerical Results and Discussions 

In our numerical calculation we take following input parameters [271 [25]: 

mw = 80.403 GeV, mz = 91.1876 GeV, 

mt = 172.5 GeV, nib = 4.7 GeV. (4.1) 

aew{0)~^ = 137.0359991, as{ml) = 0.1176, 

The QCD renormalization scale is taken to be /i = ^o(= + ^b) if there is no other 
statement, and the running strong coupling as(^^) is evaluated at the three-loop level {MS 
scheme) with the five active flavors [27] . For the numerical calculations of one-loop integrals, we 
use LoopTools2.1[T2] package to deal with 2-, 3- and 4-point integrals. The implementations 
of the scalar and tensor 5-point integrals are done exactly by using our Fortran programs as 
used in previous works [501 131] with the approach presented in Ref.t32j. And the 6-point scalar 
and tensor integrals are evaluated by using our created programs with the expressions given 
in Refs.fH] and [15]. In Ref.|16j. T. Binoth, et al., derived an analytic expression for the 
scalar hexagon function, which is convenient for the subsequent numerical integration. We 
checked also the numerical results of 6-point scalar integrals by using two methods presented in 
Refs.[ni I15j and [TB], and confirmed the correctness of our 6-point scalar integral program. For 
example, with the Set(I) input parameters in Ref.|16j we get exactly the same numerical results 
as (1.3526 x 10^^ _,_ 4.06O8 x 10"^^ i) by using both two methods. 

During our numerical calculation, we have studied the independence of the total cross section 
involving the NLO QCD corrections of process 77 — > ttbb on the soft cutoff Ss{= 2 AEi/^/s). 
To show that independence, we depict the cross section parts (T4(= atree + '^soft + '^virtual)^ C5(= 



f^hard) and NLO QCD corrected cross section a^LO versus 5s in Figsi3]^a-b) with ^/s = 800 GeV. 
As shown in these two figures, both o"4 and 0-5 obviously depend on the soft cutoff 6s, but (Tnlo 
is independent of the soft cutoff value with the best fit average value being 7.978 fh and the 
errors being less than 1.3% and 0.5% in the 5s regions of [10~^, 10^^] and [10~^, 5 x 10^^], 
respectively. In further numerical calculation we fix 5s = 10^'^. 
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Figure 3: (a) The cross section parts 0-4 (= atree + o^soft + (^virtual), <^5(= (^hard) and the NLO 
QCD corrected cross section(fTArio) of the 77 ttbb process as the functions of the soft cutoff 
5s{= 2 /S.Et / ^/s) with n = hq = mt + rrih and ^/s = 800 GeV. (b) The enlarged plot of Figl3]J^a) 
for the NLO QCD corrected cross section((TArLo) with integration error versus 5s- 

We define the K-factor as the ratio of the NLO QCD corrected cross section and the LO 
cross sectionfil' = The numerical results of the cross section and the K-factor for the 

process 77 — > tibh are plotted in FigUU^a) and FigUl^b) respectively, with the 77 colliding energy 
^/s running from 400 GeV to 2 TeV . In Figl4l[a), the full and the dashed curves correspond to 
the LO and NLO QCD corrected cross sections separately. As indicated in Fig|4l^a), the cross 
section increases quickly as the 77 colliding energy running from 400 GeV to 2 TeV , and the 
NLO QCD correction obviously enhances the tree-level cross section in the plotted range of \/s. 
The NLO QCD corrected cross section of process 77 — > tihb with ^/s = 2 TeV can reach the 
value of 15.39 fb. Figl2||b) shows the corresponding K-factor varies from 1.70 to 1.14 as the 
c.m.s energy ^/s running from 400 GeV to 2 TeV . The analysis of the contribution parts of 
NLO QCD correction shows that in the small energy region the K-factor is enhanced due to a 
Coulomb singularity effect on the contribution from the diagrams with virtual gluon exchange 
between heavy quarks. The cross sections, atree and onlo-, and K-factors at some typical -y/i 
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Figure 4: (a) The LO and NLO QCD corrected cross sections for the process 77 — > ttbb as the 
functions of c.m.s. cohiding energy(y^), (b) the corresponding K-factor versus ^/s. 





atree{fb){GRACE) 


Cftreeifb) 


CTNLoifb) 


K-factor 


500 


1.4453(4) 


1.4458(5) 


2.24(1) 


1.55 


800 


5.991(3) 


5.990(4) 


7.96(4) 


1.33 


1000 


8.001(6) 


8.000(7) 


10.18(6) 


1.27 


2000 


13.51(2) 


13.50(2) 


15.39(9) 


1.14 



Table 2: The LO and NLO QCD corrected cross sections, K-factors for 77 — > ttbb process 
with ^ = 500 GeV, 800 GeV, 1000 GeV, 2000 GeV, respectively. The LO cross section 
values obtained by using GRACE2.2.0 system are listed there for comparison. The number 
of the Monte Carlo events is 5 x 10^. 



points can be read out from FigslD[|a-b) and are listed in Table [2j For the correctness check 
of the calculation of the LO cross section for 77 ttbb process, the results obtained by using 
GRACE2.2.0 system are presented there too. 

The renormalization scale dependence of both the LO and NLO QCD corrected total cross 
sections for the process 77 — > ttbb with y/s = 800 GeV, is plotted in FigjSl In this figure the 
scale is parameterized as n/fJ-o (/^o = n^t + rnb) and the full-line and dashed-line correspond to 
the LO and the NLO QCD corrected cross sections, respectively. We can see from the figure that 
in the region 0.75 < h/hq < 4 the NLO QCD correction obviously improves the independence 
of the renormalization scale Therefore, we can conclude that the uncertainty of the cross 
section for process 77 — > ttbb due to the variation of renormalization scale /i, can be reduced by 
considering the NLO QCD corrections. 

The distributions of the transverse momenta of top- and bottom-quark(p^ and p^) with the 
colliding energy ^/s = 800 GeV, are depicted in Figl6]J^a) and Figl6]J^b) separately. In Figl6]^a), we 
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Figure 5: The LO and NLO QCD corrected cross sections for the 77 
functions of renormahzation scale /i//io(Mo = "^-t + "T-fo = 147.2 GeV). 



ttbb process as the 



can see that the NLO QCD correction obviously enhances the differential cross section predicted 
in the SM at tree- level when the pip value is less than 250 GeV. Our analysis shows that in the 
region 20 GeV < < 250 GeV, the large correction to the distribution of comes mainly from 
the contribution of the hard gluon emission process. But when pip > 250 GeV, the correction to 
the distribution of pip from the hard gluon emission process, is largely cancelled by the negative 
correction from virtual gluon and soft gluon emission contributions, then the pip distribution 
corrections become to be much smaller. That means the absolute correction (I '^'^J^t^ ~ ^^^1) 
is enhanced at low pip and reduced at large pip due to the momentum balance between top 
particles and gluons radiated from top quarks at the NLO, which reduces the momenta of the 
top quarks. For b-quark, we can see from Figl6l|b) that the large enhancement of the QCD 
corrected differential cross section(^^P), which can nearly double the LO differential cross 
section somewhere, can be appeared in the p^ value range between 20 GeV and 120 GeV, while 
in the range of pp > 150 GeV the NLO QCD correction becomes to be very small. Similar with 
the discussion for the distribution of pip, that is also the consequence of the momentum balance 
between bottom particles and gluons radiated from bottom quarks at the NLO. 



V. Summary 

In this paper we calculate the complete one-loop QCD corrections to the process 77 ttbb at 
a photon-photon collider. We present the dependence of the NLO QCD correction of process 




p^(GeV) pJ(GeV) 



Figure 6: The distributions of the transverse momenta of top- and bottom-quark for the 77 — > 
tibb process with ^/s = 800 GeV. (a) for top-quark, (b) for bottom-quark. 

77 — > ttbb on colhding energy ^/s in the SM, and find that NLO QCD correction can generally 
increase the LO cross section. It shows that in the 77 colliding energy range of 400 GeV < 
^/s < 2 TeV, the corresponding K-factor goes down from 1.70 to 1.14. We find that the 
NLO QCD correction can obviously improve the independence of the cross section for process 
77 — > ttbb on the renormalization scale fi, and the NLO QCD correction also changes obviously 
the distributions of transverse momenta of the final top- and bottom-quark states. 

Acknowledgments: We would like to acknowledge professor C.-S. Li for bringing our atten- 
tion to this issue. This work was supported in part by the National Natural Science Foundation 
of China, the Education Ministry of China and a special fund sponsored by Chinese Academy 
of Sciences. 



References 

[1] Tevatron Electroweak Working Group(for the CDF and DO Collaborations), Fermilab-TM- 
2347-E, TEVEWWG/top 2006/01, CDF-8162, DO-5064, |hep-ex/0603039/ l. 

[2] I. Bigi, H. Krasemann, Z. Phys. C7 (1981) 127; J. Kiihn, Acta.-Phys.-Austr. (Suppl.) 
XXIV (1982) 203; I. Bigi, Y. Dokshitzer, V. Khoze, J. Kiihn and P. Zerwas Phys. Lett. 
B181 (1986) 157. 

[3] C. Farreh, A. H. Hoang, Phys. Rev. D74 (2006) 014008. 



[4] V. Fadin, V. Khose and T. Sjostrand, Z. Phys. C48 (1990) 613. 

[5] M. Strassler and M.E. Peskin, Phys. Rev. D43 (1991) 1500. 

[6] W. Kwong, Phys. Rev. D43 (1991) 1488. 

[7] Y. Sumino, K. Fujii, K. Hagiwara, H. Murayama, C.K. Ng, Phys. Rev. D47 (1993) 56. 

[8] F.J. Yndurain, Phys. Lett B321 (1994) 400. 



[9] Lah Chatterjee, Cheuk-Yin Wong, |hep-ph/9501218 
[10] Han Liang, Ma Wen-Gan, Yu Zeng-hui, Phys. Rev. D56 (1997) 265. 
[11] A. Denner, S. Dittmaier, M. Strobel, Phys. Rev. D53 (1996) 44. 

[12] Chen Hui, Ma Wen-Gan, Zhang Ren- You, Zhou Pei-Jun, Hou Hong-Sheng, Sun Yan-Bin, 
Nucl. Phys. B683 (2004) 196. 

[13] A. Aepph, G.J. van Oldenborgh and D. Wyler, Nucl. Phys. B428(1994) 126, 



hep-ph/9312212| W. Beenakker, F.A. Berends, A.P. Chapovsky, Nucl.Phys. B548(1999)3. 

[14] A. Denner and S. Dittmaier, Nucl. Phys. B734 (2006) 62. 

[15] A. Denner, Fortschr. Phys. 41(1993)4, 307. 

[16] T. Binoth, G. Heinrich and N. Kauer, Nucl. Phys. B654 (2003)277. 

[17] A. Denner, S. Dittmaier, M. Roth, L.-H. Wieders, Phys. Lett. B612 (2005) 223. 

[18] T. Hahn, Comput. Phys. Commun. 140 (2001)418. 

[19] T. Hahn, M. Perez- Victoria, Comput. Phys. Commun. 118 (1999)153. 

[20] E. Byckling and K. Kajantie, 'Particle Kinematics', John Wiley and Sons, London, 1973; 
Stefan Weinzierl, |hep-ph/0006269l 



[21] J. Fujimoto, et al., Comput. Phys. Commu. 153(2003)106; T. Ishikawa, Y. Kawabata, et 
al., Minami-Tateya collaboration, 'GRACE User's Manual', Version 2.0. 



[22] B. W. Harris and J.F. Owens, Phys. Rev. D65 (2002) 094032, hep-ph/0102128 



[23] W. Beenakker, S. Dittmaier, M. Kramer, B. Plumper, M. Spira and P.M. Zerwas, Phys. 
Rev. Lett. 87 (2001) 201805; Nucl. Phys. B653(2003)151-203. 

[24] S. Catani and M.H. Seymour, Phys. Lett. B378 (1996) 287 and Nucl. Phys. B485 (1997) 
291 [Erratum-ibid. B 510(1997) 291]. 

[25] S. Catani, S. Dittmaier, M.H. Seymour and Z. Trocsanyi, Nucl. Phys. B627 (2002) 189. 

[26] L. Reina, S. Dawson and D. Wackeroth, Phys. Rev. D65(2002)053017; P. Wu, W.-G. Ma, 
R.-Y. Zhang, Y. Jiang, L. Han, and L. Guo, Phys. Rev. D73, (2006)015012. 

[27] W.M. Yao, et al., J. of Phys. G33,l (2006). 



[28] F. Jegerlehner, DESY 01-029, hep-ph/0105283 



[29] E. Boos, V. Bunichev, et al., (the CompHEP collaboration), Nucl. lustrum. Meth. A534 



(2004) 250-259, 'hep-ph/0403Tl3 



[30] Y. You, W.-G. Ma, H. Chen, R.-Y. Zhang, Y.-B. Sun, H.-S. Hou, Phys.Lett. B571 (2003) 
85-91. 

[31] R.-Y. Zhang, W.-G. Ma, H. Chen, Y.-B. Sun, H.-S. Hou, Phys.Lett. B578(2004) 349-358. 
[32] A. Denner and S. Dittmaier, Nucl. Phys. B658(2003) 175. 



